Flocking dynamics and mean-field limit 
in the Cucker-Smale-type model 
with topological interactions 

Jan Haskovec^ 

' Abstract. We introduce a Cucker-Smale-type model for flocking, where the strength of interaction be- 
. tween agents depends on their topological distance, measured in units of agents' separation (number of 
CN ! intermediate individuals separating the agents), which makes the model scale- free. This is motivated by 
■ recent extensive observations of starling flocks, which suggested that interaction ruling animal collective 
behavior depends on topological rather than metric distance. We study the conditions leading to asymp- 
totic flocking in the topological model, defined as the convergence of the agents' velocities to a common 
vector. The shift from metric to topological interactions requires development of new analytical methods, 
taking into account the graph-theoretical nature of the problem. Moreover, we offer a purely formal 
derivation of the mean-field limit of large populations, recovering kinetic and hydrodynamic descriptions. 
In particular, we introduce the novel concept of topological distance in continuum descriptions, which 
is applicable to a broad variety of models of collective behavior. As an example, we shortly discuss a 
topological modification of the attraction-repulsion model and illustrate with numerical simulations that 
the modified model produces interesting new pattern dynamics. 
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O ■ 1 Introduction 
^ ! 

■ Collective behavior of large animal groups with self-organization into robust complex patterns is a fas- 

, cinating natural phenomenon [19, 16]. Prominent examples are schools of fish, which can move in a 
• • I rather ordered formation, change direction abruptly, or, under predator threat, swirl like a vehemently 

• stirred fluid [20]. Flocks of hundreds of starlings can fly as a uniformly moving group, as well as produce 
^ . turbulent aerial displays [11]. Apart from its biological and evolutionary relevance, collective phenomena 

I play a prominent role in many other scientific disciplines, such as robotics, control theory, economics and 
social sciences, see, e.g., the recent surveys [7, 22]. 

Regarding its pratical and theoretical importance, it is not surprising that the topic attracted wide 
attention of physical and mathematical scientific communities. Many different types of mathematical 
models have been proposed and studied during past decades. Their common assumption is that per- 
manently moving self-propelled individuals (agents) interact and adapt their behavior according to their 
conspecifics. Among them, the Cucker-Smale model is very well known, in particular due to its relative 
simplicity and the implied convenience for analytical study. In fact, the original motivation in [8, 9], 
where the model was introduced, was to describe language evolution. Only in subsequent studies it has 
been eventually related mainly to the description of the emergence of flocking in groups of interacting 
animals. The model considers a group of G N agents with time-dependent locations Xi{t) G M.'^ and 
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velocities Vi{t) G W^, subject to 



Xi = Vi, (1.1) 

^ jV X] ~ ' for i = 1, . . . , A^. (1.2) 

In the classical setting, the so-called communication rates gij depend on the metric distance between the 
z-th and j-th. individual, gij = g(\xi — Xj\) with a fixed function g : [0,oo) — )■ [0, oo). In particular, the 
rate function considered in [8, 9] and most of the subsequent papers is 

= (^2 + ,2)/3 > (1-3) 

with A, a and P positive parameters. Then, it was shown that if /3 < 1/2, the model exhibits the so-called 
unconditional flocking, where for every initial configuration the velocities Vi(t) converge to the common 
consensus value Xli^i Vi{0) as t — )■ oo. On the other hand, with /3 > 1/2 the flocking is conditional, i.e., 
the asymptotic behavior of the system depends on the values of A and a and on the initial configuration. 
This result was first proved in [8, 9] using tools from graph theory (spectral properties of graph Laplacian), 
and slightly later reproved in [21] by means of elementary calculus. Another proof has been provided in 
[15], based on bounding (1.1)-(1.3) by a system of dissipative differential inequalities, and, finally, the 
proof of [6] is based on bounding the maximal velocity. Moreover, the case of singular communication 
rate g{s) = X/s^ was studied in [15]. 

Another classical topic of study is the derivation and analysis of the mean-field limit as N ^ oo, 
which leads to the Vlasov-type kinetic equation 

dt + V^- {vf) + V„ • {g[f]f) = 0, (1.4) 

with 



0[f]it,x,v) = - / g{\x -y\){v -w)f{t,y,w)dwdy, 



where / = f{t,x,v) is the agent distribution function; see, e.g., [4, 6, 15]. A stochastic version of the 
model was considered in [14]. 

In our paper we introduce a modification of (1.1)-(1.2) where the communication rates gij depend 
on the topological, rather than metric, distance between the agents i and j. This is inspired by recent 
extensive observations of starling flocks in [2] , which indicate that starling interactions are scale-free and 
their strength depends on the topological distance between individuals, measured in units of average bird 
separation. In other words, the relevant quantity is how many intermediate individuals separate two 
birds, not how far apart they are in the metric sense. Moreover, it is the very shape of the interaction 
that depends on the topological distance, not simply the cut-off or the range. Consequently, we introduce 
the (normalized) topological distance between the i-ih and j-ih. agent as 



N 



(^ij = < k<N-\xk-xj\ < \xi-xj\} = j^^xi^^-::^) , (1-5) 

k=l 



"J I 



where x is the characteristic function of the open unit ball in M*^. We will then consider the model (1.1)- 
(1.2) with the communication rates gij depending on the topological distance aij through gij = g{aij), 
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where g : [0,1) [0, oo) is a nonnegative function (note that < aij < 1). In our paper, we are 
interested in studying the asymptotic flocking behavior of this model in dependence of the function g 
and the configuration of the agents. We will show that fundamental for the topological model to exhibit 
flocking asymptotic is the connectivity of the graphs of inter-agent interactions. Moreover, we will formally 
derive the mean-field limit of the model, which will be of the Vlasov-type (1.4). Finally, we will propose 
that the idea of introducing scale-free interactions is applicable to a broad spectrum of discrete and 
continuum models of collective behavior, and mention the well-known attraction-repulsion model [10] as 
a particular example. 

Before we proceed to the mathematical analysis of the topological Cucker-Smale model, let us make a 
short remark about the biological relevance of topological interactions, as opposed to metric ones. Indeed, 
as argued in [2], the difference between metric and topological interactions has a major impact on the 
global flocking behavior of the corresponding models, and, consequently, on their biological relevance. The 
most important evolutionary advantage of collective behavior is avoidance of predation, for which strong 
cohesion of the group is essential [20]. Numerical experiments performed in [2] and [3] strongly suggest 
that only topological interactions grant such robust cohesion and, therefore, high biological fitness. The 
intuitive explanation of this fact is the high "flexibility" of topological interactions, compared to metric 
ones: In the metric case, interaction effectively vanishes and cohesion is lost whenever the interindividual 
distance becomes larger then the prescribed radius. In contrast, topological interaction stays effective also 
over long distances or in the case of low individual density, and so keeps the cohesion even in the presence 
of strong perturbations, of which predation is the most relevant. Another work revealing the qualitative 
difference between the two kinds of interaction is [12], where a topological Voronoi Vicsek-type model 
exhibited an ordered moving phase with novel long-range correlations, even though no direct long-range 
connections emerged, and the transition to collective motion exhibited critical properties different from 
known universality classes. 

The obvious biological relevance of topological interactions contrasts with the fact that the mathe- 
matical flocking community so far focused almost exclusively on models with metric interactions. This is 
probably best explained by their accessibility to mathematical analysis and studies of global quantitative 
properties. In our paper wc aim to make a step towards closing this gap by studying the flocking proper- 
ties of the discrete topological Cucker-Smale model and deriving its kinetic and hydrodynamic limits. The 
paper is organized as follows: In Section 2 we offer a detailed formulation of the model, relate it to known 
models of control theory and discuss its main properties. We also provide two examples pointing out that 
the connectivity of the graphs of inter-agent interactions play a central role for the asymptotic flocking 
behavior. In Section 3 we provide our main result about the unconditional flocking. In Section 4 we 
derive the formal kinetic and hydrodynamic limits of the model as the population size N tends to infinity 
and briefly discuss their flocking properties. Finally, in Section 5 wc argue that our approach is applicable 
to a broad class of flocking models. As a particular example, we mention the attraction-repulsion model 
[10] and illustrate with numerical simulations that the modified model produces interesting new pattern 
dynamics. 
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2 Formulation of the model and its main properties 



We consider the Cucker-Smale-type model 

Xi = Vi, (2.1) 
1 ^ 

Vi = — y29ioiij)ivj-Vi), for 1 = 1,..., N, (2.2) 

with the topological distance aij defined by (1.5). Clearly, the word "distance" is a misnomer, since the 
mathematical definition of a distance is violated in two ways: the symmetry and triangle inequality do 
not hold. However, we will allow ourselves this abuse of terminology in this paper. At least, the definition 
(1.5) guarantees that aij = if and only if Xj = Xj. 

The behavior of the system is determined by the rate function g : [0, 1] — > [0, oo). Of course, for a 
fixed N, we only need the values of g on the set {0, 1/A^, . . . , {N — 1)/A^}. However, in order to being 
able to formulate results independently of N and later perform the limit as — t- oo, we will assume 
that g is defined on the whole interval [0,1]. This will also allow for using a smoothened version of 
the characteristic function x in (1-5) if necessary. Prom the modeling point of view, it makes sense to 
assume that interactions with closer (in the sense of the topological distance) agents are stronger, which 
corresponds to an nonincreasing g; however, we do not require the validity of this assumption for the 
forthcoming analysis. 

Note that the right-hand side in (2.2) is scaled by instead of the scaling by 1/N in (1.2), with 

N 

7iv = J]5((^-l)/A^)■ (2.3) 

i=l 

This is motivated by the fact that Ylf=i di^^ij) = In for every i = 1, . . . ,N, so that ^ Z^jLi 9{(^ij)vj is 
a convex combination of the vectors vj. The same type of scaling was considered in [17], going under the 
name of relative distance; however, it is still a metric distance, with gij = g{\xi — Xj\). The advantage 
of the model [17] is that it does not involve any explicit dependence on the number of agents; just 
their geometry in phase space is taken into account. Therefore, it cures the drawback of the classical 
Cucker-Smale model, where the motion of an agent is modified by the total number of agents even if its 
dynamics is only influenced by essentially a few close neighbors. Our approach, which is to replace the 
metric distance by its topological counterpart, treats this drawback in a much more radical way, and, as 
explained in Section 1, has the advantage of being supported by actual observations of biological collective 
motion (starlings). 

The mathematical disadvantage of introducing relative or topological distances is that the model 
looses its symmetry. Indeed, observe that in the classical Cucker-Smale model (1.1)-(1.2) we have gij = 
g{\xi — Xj\) = gji. Symmetry is the cornerstone for studying the long time behavior of its solutions, since 
it implies that the total momentum V{t) := jj Vi{t) is conserved, and with appropriate assumptions 
one can prove the decay of the fluctuations of the velocities about V{t) = V{0). In fact, the necessary 
condition for momentum conservation is that the communication matrix G = (gij) G M^^^ is balanced, 
i.e., its row and column sums are equal, 

N N 

^9ij = ^9ij- 

i=i j=i 
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Neither the model of [17] nor our topological model is balanced. Therefore, new analytical techniques to 
study the flocking behavior need to be developed. The approach of [17], based on the notion of active 
sets, requires the rather strong assumption that all the communication rates gij are strictly positive, i.e., 
there is a direct communication between every two agents. In our approach, we admit the more realistic 
setting where some or many of gij actually vanish, and we will show that the strong connectivity of the 
underlying directed graph is sufficient and also necessary for asymptotic flocking to take place. 

For a fixed N, let us rescale the model such that 7jv = 1. Then (2.2) can be written in the form 

v = {G-I)v = -Lv , 

where the rows of G M^^'^ are composed of the vectors vi, . . . ,vn, and G = (gij) € M^^^ with 
gij = g{oiij) is a stochastic matrix (in the sense that all row sums are equal to 1). The matrix L := I — G 
is then the discrete Laplacian corresponding to the directed graph (called digraph in the sequel) spanned by 
the agents with edge weights given by gij > 0. In this sense we will indetify the topological configuration 
of the system with the matrix G or, equivalcntly, L. Obviously, with N agents one can have at most 
A^! different topologies G'' , k = 1,...,N\, and the system switches between (some of) them during its 
temporal evolution. Because the trajectories Xi of the agents are continuous, in any time interval there 
can be only countably many switches and any two consecutive switches are separated by a nonempty open 
time interval, where the configuration does not change. We can then understand the system (2.1)-(2.2) 
as a switching communication network, see, e.g., [18], 

x{t) = v{t), (2.4) 
v{t) = -L'^(^W)i;(t), (2.5) 

with -L'^^^^*)) = - / and the function a : M^'^'' ^ {1, • • ■,N\} switches between the topologies 

according to the current configuration x{t). We introduce the following notion of solution: 

Definition 1 Let the curve {x{t),v{t))t>o G M.^^'^ x M.^^^ be globally continuous. Denote a : M.^^'^ — > 
{1, . . . , A^!} a switching function with at most countably many switches on (0, oo), and denote {Ik}kef^ the 
system of open intervals such that a = const, on every and (JfeeN-^fc ~ °°)- '^'^^^ continuous 
curve {x{t),v{t))t>Q a solution to (2.4)-(2.5) if it solves the differential equation on every open interval 
Ik, keN. 

For the case of a fixed configuration G = (gij) with a strongly connected digraph (let us recall that 
strong connectivity means that every pair of vertices u, w is connected by a directed path from u to w 
and a directed path from it; to u). Corollaries 1 and 2 of [18] imply that an asymptotic consensus is found 
even if G is not balanced: 

Proposition 1 Consider the communication network v = {G — I)v with a fixed topology G that is a 
strongly connected digraph. Let ^ = (^i, . . . ,CAf)^ be a nonnegative left eigenvector corresponding to the 
zero eigenvalue of L = —{G — /), i.e., $,^L = 0, such that Yl^i^i > 0- Then, an asymptotic consensus 
v°° G is reached exponentially fast as t oo, i.e., 

lim Vi{t) = v°°, for 1 = 1,..., N, 
with v°° being the convex combination of Vi, 



The proof [18] is based on the simple observation that ^"^v is an invariant quantity (note that the graph 
Laplacian has always a zero eigenvalue). Let us point out the two fundamental assumptions of the above 
Theorem: (i) the topology is fixed and (ii) the corresponding digraph is strongly connected. Of course, the 
above result would be applicable for our system (2.4)-(2.5) as soon as we knew that the system undergoes 
only a finite number of switches. We would then simply pick the last attained configuration and if it 
was strongly connected, wc would conclude that an asymptotic consensus will be reached exponentially 
fast as t — >■ oo. However, the following example shows that a special choice of the communication rate 
g together with a particular initial configuration leads to an infinite number of switches. Even worse, it 
may be possible that the system (2.1)-(2.2) does not admit global solutions: 

Example 1 Let us consider a group of 7 agents, denoted by i = —3, —2, —1, 0, 1, 2, 3, moving on the real 
line, with initial positions 

x_i(0) = -6, x_2(0) = -9, x_3(0) = -10, 
xi(0) = 6, X2(0) = 9, X3(0) = 10, 

and initial velocities Vi = —1 for i = —1,-2, —3 and Vi = 1 for i = 1,2,3. For the agent i = we 
prescribe the initial datum 

xo(0) = 0, vo{0) = c, (2.7) 

for some c > 0. Let us consider the model (2.1) -(2.2) with weight g given by 

5(2/7) = 1, gii/7)=0 fori ^2, 

i.e., evert agents interacts exclusively with its second closest neighbor. Let us assume for a moment that 
we are able to choose c > such that —1 < XQ^t) < 1 for all t > 0; we will justify this assumption later. 
Then, the agents from the two triplets i € {—1,-2,— 3} and i G {1,2,3} will interact only with agents 
from the same group, so that the velocities of these six will remain constant and equal to their initial 
values. The agent i = will interact with i = —1 if Xo{t) > and with i = 1 if xo{t) < 0. Therefore, 

xo = vq, vo = -sign(a;o) - vq. (2.8) 

Due to the initial condition (2.7), the agent first moves to the right, so that sign(xo) = 1. Then (2.8) is 
resolved by 

xo(t) = (c+l)(l-e-*)-t, ^;o(t) = (c + l)e-* - 1, (2.9) 

Note that since Xg(0) = c > and linit^^xo{t) = —oo, there exists a positive time t = r(c) such 
that xq{t(c)) = 0, i.e., the agent will return to the origin eventually, and then its trajectory will be be 
subject to (2.8) with sign(xo) = —1. The "turning point" tt, when vo{tt) = 0, is tt = ln(c + 1) and 

Xo{tt) = c — ln(c + 1), so for c small enough the agent will be confined to the strip —1 < a;o < 1, which 
verifies the above assumption. The connectivity diagram of the two configurations with xq > and xq < 
is visualised in Figure 1. 

However, the agent will never exhibit a periodic trajectory, because vq{t{c')) > —c. Indeed, since 
xait) + voit) = c — t, we have vq(t{c)) = c — r(c). Now if there was a c> such that vo{t{c)) = —c, this 
would imply r(c) = 2c and inserting this into (2.9), we would get 

xo{t) = (l + r/2)(l-e-^)-T, 
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Figure 1: Connectivity diagram for Example 1, upper panel for xq > and lower panel for xq < 0. An 
arrow is pointing from Xi to Xj if and only if gij > 0. 

which, as can be easily checked, is equal to zero only for r = 0. Therefore, the movement of xq is never 
periodic, but still the agent crosses the origin x = infinitely many times, thus changing the configuration 
of the system infinitely many times. To see this, let us calculate the implicit derivative of xo{t{c)) = 
with respect to c, 

T (c) 



(c + l)e-'^W - 1 c - r(c) 

Realizing that vo{t{c)) = c — r(c) < (the agent can only return back to the origin with a nonpositive 
velocity), we find t'{c) > for all c > 0. Moreover, using the L'Hospital rule, we calculate 

r'(0) = 2, r"(0) = -4/3. 

This implies that for c small enough, vo{t(c)) = c—t{c) > —c, i.e., the speed s{c) := |fo(r(c))| of the agent 
when returning to the origin is smaller than its initial speed c when it was leaving the origin. Moreover, 
the time needed to return to the origin, t{c), satisfies c < r(c) < 2c for all c > 0. A plot of the trajectories 
and velocities of the agents is provided in Fig. 2. 



Thus, the agent follows a nonperiodic trajectory, starting from xo(0) = and vo{0) = c > 0, and, 

since s{c) > for all c > 0, returning to the origin infinitely many times, at instances t(c), r(s(c)), 
t{s(s(c))) , . . . , with speeds, respectively, s(c), s(s(c)), s{s{s{c))), .... Consequently, the configuration of 
the system switches infinitely many times. However, the fact that both the sequence of return times and 
the sequence of speeds are strictly decreasing ( and tend to zero, since otherwise a periodic solution would 
exist) immediately poses the question whether the series of return times diverges or not. If it does diverge, 
then the nonperiodic solution of the ODE (2.8) exists globally in time. On the other hand, a convergent 
series of return times leads to a local solution only. 

Unfortunately, we were not able to find an analytic answer to this question; let us remark that since 

lim ^ = lim = r'(0) - 1 = 1, 

C c-*-0 C 

the ratio test for convergence of series is inconclusive, and we were not able to obtain the answer from any 
of the higher- order tests. The only trivial observation is that the answer does not depend on the particular 
value of c. 




Figure 2: Agent locations (left panel) and velocities (right panel) in Example 1. 

Our second example shows that even though only a finite number of switches takes place, strong 
connectivity of the digraph of the last attained configuration is necessary for a consensus to be found. 
Indeed, if the digraph is only weakly connected (let us recall that a digraph is called weakly connected 
if replacing all of its directed edges with undirected edges produces a connected undirected graph), the 
system may not find any consensus: 

Example 2 Let us consider a group of 7 agents, denoted by i = —3, —2, —1, 0, 1, 2, 3, moving on the real 
line, with initial positions 

x_i(0) = -1, a;_2(0) = -3, x_3(0) = -6, 
xi(0) = l, X2(0) = 3, X3(0)=6, 

and initial velocities Vi = —1 for i = —1,-2,— 3 and Vi = 1 for i = 1,2,3. For the agent i = we 
prescribe the initial datum 

xo(0) = 0, vo{0) = 0. 

Let us consider the model (2.1) -(2.2) with weight g given by 

5(1/7) =5(2/7) = 1, 5(V7)=0 /ori>3, 

i.e., every agents interacts with its two closest neighbors, updating its velocity according to the arithmetic 
mean of their velocities. 

Instead of attempting to find an analytic solution to the corresponding ODE system, we performed a 
numerical simulation (see Fig. 3), showing that the system takes two configurations: 

• First, for t < to, with to ~ 10, the system stays in the initial configuration with a strongly connected 
digraph, plotted in the upper panel of Fig. 4- The configuration remains symmetric with xq = 0, 
vo = ^. 
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Figure 3: Agent locations (left panel) and velocities (right panel) in Example 2. 



• After time to, the distance between xi and becomes smaller than the distace between xi and xq. 
Due to the symmetry, the same happens for x-i and x-^. Therefore, the connectivity of the digraph 
changes from strong to weak, see lower panel of Fig. 4- Then, the agents of the group {1,2,3} 
interact only among themselves, and the same holds for the group {—1, —2, —3}. Therefore, each 
of the two groups will find their own velocity consensus as t ^ oo, which will be nonzero velocities 
with opposite signs. The zeroth agent will stay with xq = 0, vq = 0. Consequently, no global velocity 
consensus will be achieved. 

Let us note that in this example the initial digraph was strongly connected, but this property got lost during 
the evolution (due to switching to another configuration). 



3 Asymptotic flocking 

We now present our main result, which applies to the full model (2.1)-(2.2), where we assume global 
existence of solutions. In fact, in the proof we do not explicitly make use of the fact that the function 
a{x) switches the topological configurations according to the topological distances gij = g{aij). The 
only relevant property is that there is at most a countable number of switches, which is guaranteed by 
the continuity of the trajectories x{t). The fundamental assumption is then that there exists a topological 
configuration with a strongly connected digraph where the system spends an infinite amount of time. This 
of course is a strong assumption and might be difficult to verify practically, however, the above Example 
2 shows that it is necessary. Prom this point of view, the forthcoming Theorem is optimal. 

Theorem 1 Assume that the system (2.1)^(2.2) has a global solution {x{t),v{t))t>o. Moreover, assume 
that there exists a topological configuration, say = {g^j), with a strongly connected digraph, where the 
system spends an infinite amount if time, i.e., 

|{t > 0; g{aij) = g\ for all i,f = \, . . . ,N}\ = +cx). 
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Figure 4: Connectivity diagram for Example 2. Upper panel shows the strongly connected configuration 
until t ~ 10, which then changes to the weakly connected configuration in the lower panel. An arrow is 
pointing from Xi to Xj if and only if gij > 0. 

Then the system finds an asymptotic velocity consensus, i.e., there exists a vector v°° G such that 

lim ='u°°, for alii = 1,...,N . (3.1) 

t-^oo 

Proof: The proof will be carried out in several steps. 

Step 1: Convergence of maximal speed. Refering to Definition 1, let us denote by 1^ := {tk-i,tk), 
k £ N, the system of disjoint open time intervals where the system does not change its topological 
configuration, and such that UfeeN-^fc ~ [0, oo). Inspired by [6], let us define the function 

ujit) := max \vi(t)\ for alH > , 
^ ' i=i,...,N ' ^ ^' — ' 

and, moreover, denote M{t) := argmax^^j^ Arj^i(OI; if M{t) is not uniquely determined (i.e., there are 
several vectors Vi(t) with maximal length), we choose one of the indices arbitrarily, but in such a way 
that M(t) stays constant on the longest time interval. Since the number of particles N is finite and their 
trajectories in the (x, i')-space are continuous, there exists an at most countable system of open disjoint 
intervals {]Ck)kef>i, such that UfceN^fc ~ [0, oo), and M{t) is constant on every /C^. To ease the notation, 
we will usually skip the explicit dependence of M on t (or k) in the sequel. 

By intertwining the two systems {Ik)kef>i ^^'^ (^fc)feeN5 we construct another at most countable system 
of disjoint intervals, denoted by abuse of notation again by {Ik)keN, such that the topological configuration 
does not change and the index M is constant on each Ik- Then, on every Ik the trajectories of the particles 
in phase space are smooth and we can write 

^"^'^^'^ " ^3Mj{Vj - Vm) ■ VM 
N 

< y^,9Mj {\vj\ - \vm\) \vm\ ■ 
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Dividing by \vm\ (note that if \vm\ was zero, there would be nothing to prove), we obtain 
d d ^ 

-^t^(i) = -^\'"m\ < ^9Mj i\vj\ - \vm\) < on every Ik, 

where the nonpositivity of the right-hand side is due to the maximahty of \vm\- Consequently, uj{t) is a 
globally continuous, nonincreasing and nonnegative function, so that there exists an < Wqo < ^{^) such 
that limt^^Lu{t) = Uoo- 

Step 2: Convergence of one velocity. Since, by assumption, the system spends an infinite 
amount of time in the strongly connected configuration we can pick the corresponding subsystem 

out of (Xfc)jfcgpj, with infinite length. By a further subselection we get the system {Tkn)n&i of infinite 
length, where M{t) = Mq for some fixed 1 < Mq < A^. Therefore, denoting := UneN-^^n' have the 
configuration and the maximal vector index Mq for all t ^I^. Moreover, 

\vMo{t)\ ^ i^oo and -^\vMo{t)f^O as t ^ cx), t G 1°, 

where the convergence of the time derivative is due to the monotonicity of {vmo (t) \ ■ Since the configuration 
G° is fixed and strongly connected, there exists an index jo such that QmoJo ='■ > on iP . Then we 
have, for t G X°, 

Id ^ 

where we used the inequality {vj — vmq) ' vmq < implied by the maximality of I^MoI- Now, since the 
left-hand side tends to zero as t ^ oo and is a constant, we have 

'"jo • VMoit) - \vMo{t)f -^0 as i ^ oo, t G X°, 

and this futher implies vjq ■ VMoit) — ^ f^^- Finally, since vjq ■ vmq = and only if vj^ = vmq (equality 

in the Cauchy-Schwartz inequality), we have 



(vjo - VMo){t) -^0 ast^oo, tel^. (3.2) 



Moreover, we calculate 



Id ^ 

2d^l^.7o(0P = ^9jo,li'"l-'"3o)-'"joit) 
1=1 



N 



X] Sjo,l [{'"I - ^Mo) • VMo + (vi - VMo) ■ (vjo - VMo) + (vMo - Vjo) ■ Vj^] 



1=1 

N 



^ 5Z Sjo,l ivi -VMo)- VMo (t) + (0) \vjo - VMo I 



1=1 



< 3u{0)\vjQ - VMo\, 



11 



where we used the estimate maxd?;^^!, \vi\) < uj{0) and the maximahty of vmq- By (3.2) we have then 

Hmsup -T'\vjo{t)\'^ < 0. 

t^oo,t£l° 

Since li'jolOl ~^ ^oo from below on X° as i — )■ oo, wc conclude that 

where is a system of subintervals of X", still of infinite Lebesgue measure. 

Step 3: Convergence of all velocities. We will show that (3.2) holds for any index j G {1, . . . , A^}. 
Due to the simple connectivity of the digraph there exists a path Mq jo '-^ ji ^ ■ ■ ■ '-^ je 3, 

such that (/Mojo > 0> ffjoji > 0, . . . , gj^j > on 1°. We proceed inductively, showing first that the results 
(3.2) and (3.3) hold for ji as well. Indeed, passing to the limit in 

Id ^ 

2 I ^io (*) I ^ < Yl Sjo,l (vi - VMo) ■ VMoit) + 3uji0)\ - VMo \ 
1=1 

< 9 jo ,h ('"h - VMo)- VMo (*) + 3a; (0) I Vj^ -vmoI 
we obtain, due to (3.3) and the maximality of vmoi 

< lim {vj^ - VMo) ■ VMoit) < 0, 

t->oo,t6Xi 

which immediately gives (3.2) for vj^ on I^. Using this result, we argue as before to conclude 

lim -^b,-(i)|2 = o 

with a system of subintervals of of infinite Lebesgue measure. This is (3.3) for vj^ on I^. Proceeding 
inductively, after a finite number of steps we reach the index j. 

We conclude that there exists a sequence {tk)k&n C X", tk — )• oo, such that, for all j = 1, . . . ,N, 

Vj{tk) - VMoitk) ^ ask^oo. 

Step 4: Convergence of scalar products. We define the function P(t) to be the minimal scalar 
product of the velocity vectors at time t>0, 

Pit) ■■= . min Viit)-Vjit), 

i,]=l,...,N 

and denote a(t) and 6(t) the indices of the two vectors with the minimal scalar product. Note that we 
do not exclude the possibility ait) = bit). If a and b are not uniquely determined, we choose their values 

arbitrarily from the present options, but such that their values remain constant on the maximal time 
interval. As before, there exists an at most countable system of open disjoint intervals iJ'k)kef^ such that 
the topological configuration does not change and the indices ait) and 6(i) are constant on each jTfe. Then, 
for t & Jk^^^ have 

d d ^ AT 

-^^(*) = ■ = m 3ajiVj -Va)-Vb + Y^ gbjiVj - Vb) ■ Va > 0. 

3=1 j=l 
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Due to the global continuity of P{t), it is a nondecreasing function, bounded above by oj'^{t). Consequently, 
there exists a finite limit of /3{t) as t ^ oo. In particular, for the sequence tk ^ oo constructed in the 
previous step, we have 

/3(ifc) = {Va ■ Vb){tk) = {Va - VMo) ' Vb + VMq ' {Vb " VMq) + I'WMqP ^lo aS OO, 

so that linit^oo /^{t) = w^. 

Step 5: Conclusion. Since 



lim Pit) = lim u^t) = w^, 
^2 



and since f5{t) < Vj{t) ■ VM{t) <u{t) for all indices j € {1, • . • ,iV}, we obtain 

lim Vj{t)-VM{t) =ul^. 

Moreover, Vj ■ vm = only if Vj = vm, so that 

lim Vj{t) -VM{t) = 0. (3.4) 

Finally, fixing any nonzero vector ^ G M'^, we calculate 

d 



N 

^9MjiVj-VM)-^ 



AT 

<\^\^9Mj\vj{t) -VM{t)\ ^ ast^oo. 

Due to the arbitrariness of the choice of ^ and the boundedncss of \vm\, there exists a vector v°° such 
that vmq -> as t — ^ oo. This, together with (3.4), implies (3.1) and concludes the proof. 



Remark 1 The "classical" definitions of time- asymptotic flocking in the Cucker-Smale model, see, e.g., 
[8, 9, 15], pose, in addition to the velocity alignment, the requirement of uniform boundedncss of the 
position fluctuations in time (formation of a group ), 

N 

sup^^ \xi{t) — Xc{t)^ < OO, (3.5) 

where Xc{t) = '^i=iXi{t) is the centre of gravity. Such a result can be easily obtained for our topological 
model under the assumption of a finite number of switches with a strongly connected final topological con- 
figuration. Indeed, Proposition 1 provides then an exponential convergence of the velocities to a consensus 
and (3.5) follows straightforwardly. However, in the general setting with possibly infinitely many switches, 
we are only able to provide the estimate 

max \xiit)\ < max IxjfO)! + a;(0)i, 
i=l,...,N ' ^ ^' - i=l,...,iV ' ^ ^' ^ ' 

which is a trivial consequence of the bound on the maximal velocity \vm\ < a;(0). 

Remark 2 In the classical Cucker-Smale model, or any of its modifications with a balanced communi- 
cation matrix, the mean velocity V{t) = '^ii^) invariant, so that whenever asymptotic flocking 
takes place, the consensus velocity is a priori given by v°° = V{Q). In contrast to that, our model (2.1)- 
(2.2) does not seem to posses any invariants that would allow us to predict v°° from the initial datum, 
beyond the trivial fact that v°° is a convex combination of Vi{0), i = 1, . . . , N . We can therefore, similarly 
as in [17], consider v°° as an emergent property of our model, in the sense that the asym,ptotic consensus 
v°° is encoded in the dynamics of the system and not just as an invariant of its initial configuration. 
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4 Kinetic and hydrodynamic limits 



In this section we provide the formal mean-field limit of (2.1)-(2.2) as iV — )• oo. Unfortunately, we 
cannot apply the classical technique of writing down the Liouville equation for the AT-particle distribution 
function and applying the molecular chaos assumption to obtain the one-particle kinetic equation, see, 
e.g., [13]. The reason is that the interaction strength gij = g{aij) is a nonlinear function of the topological 
distance aij given by (1.5), which depends on the spatial configuration of the whole system. Therefore, 
it is not possible to obtain a closed BBGKY hierarchy and the classical method [13] fails. 

To circumvent this problem, we will formally decouple the system (2.1)-(2.2) from (1.5) by considering 

Xi — Vi , 

Vi = —y Qijivj - Vi), for i = 1, . . . , iV, 
Tat ^ 



with a prescribed (time dependent) matrix G = {gij) with all row sums X]j=i 9ij equal to 'jn- Introducing 
the empirical measure 



/^(^, ^) := ^ E S{x - Xi{t))5{v - Vi{t)) (4.1) 

i=l 

and assuming that there exists a limit — >■ / as iV — >■ oo, we obtain the mean-field Vlasov equation 
(see, e.g., [21]) 

^ + v-V.f = V,-{g[f\f), (4.2) 



with 



Q[f]{x,v) = -[ [ g{x,y){w - v)f{y,w)dwdy , 



where 7 = limjv-^^oo The undefined term g{x,y) is the formal limit of gij = g{aij) and is evaluated by 
writing the topological distance (1.5) as 



Z Xi \ fjV/ 



JwdJ^d \\Xj-Xi\J 

where % is the characteristic function of the unit ball in W^. The formal limit ss N ^ 00 gives then the 
continuum topological distance between the locations x and y in R'^, 

a[f]{x,y)= [ [ x(-^^)f{z,v)dvdz. (4.3) 
jRdJT^d \\y-x\J 

Finally, substituting g{a[f]{x,y)) for g{x,y), we obtain 

g[f]{x,v) = - / g{a[f]{x,y)){w -v)f{y,w)dwdy, (4.4) 

with 7 = /q 5((s)ds. 



14 



To carry out a rigorous derivation of the kinetic system (4.2), (4.4) in the mean-field hmit is more a 
difficult task than in the case of the classical Cucker-Smale model. In particular, due to the dependence 
of the pairwise interaction strength on the topological distance (and thus on the overall configuration of 
the system), the techniques of [21] or [4] do not carry over to our model (2.1)-(2.2). Since this topic is not 
central to this paper, we postpone the study of rigorous mean-field limit derivation and well-posedness of 
the kinetic equation (4.2), (4.4) to a future work. 

Instead, let us derive the formal hydrodynamic limit for the system (4.2), (4.4). Defining the mass 
and momentum densities 

p{t,x)= f{t,x,v)dv, pu{t,x)= f{t,x,v)vdv, 



we have 



jRd \\y-x\J 

/ G[f]{x,v)f{x,v)dv = / g{a[p\{x,y))p{x)p{y)[u{x) -u{y)]dy. 

JRd 7 JRd 

The latter expression can be further simplified by noticing that the topological distance a[p]{x,y) depends 
on y only through \x — y\, so it can be written as a[p]{x,y) = ax{\x — y\), with 

ax{r)= [ Rx{s)ds, Rx{s) = [ p{y)dS{y), 

Jo JSsix) 

where dS is the surface measure on the sphere Ss{x) with radius s > 0, centered at x. Therefore, we 
have -^axir) = Rx{r) for every fixed x eM.'^ and 

/ 9{<^[p]{x,y))p{y)dy = f g{axir))-^axir)dr = f g{s)ds = j, 
JRd Jo ur Jo 

where we used the normalization J^^ p{x) dx = ax{oo) = 1. Thus, we obtain 

Q[f\{x,v)f{x,v)dv = p{x)\u{x) -uix)] with u{x) = - [ g(a[p]{x,y))p{y)u{y)dy. (4.5) 

1 jRd 

Then, integrating the kinetic equation (4.2), (4.4) againts the moments (1, v) and taking the usual monoki- 
netic closure assumption f{t,x,v) = p{t,x)d{v — u{t,x)), we obtain the closed Euler-type system 

|^ + V-(H = 0, (4.6) 

— + {u-V)u = u-u. (4.7) 

The right-hand side in the u-equation decribes the tendency of agents with velocity u to relax to the local 
average velocity u, defined in (4.5). The word "local" is to be understood in the topological sense, i.e., 
every agent measures the average velocity of a certain portion of its closest neighbors, but independently 
of their actual metric distance. 

The Euler system (4.6)-(4.7) falls into the class studied in Section 6 of [17]. The authors adapt their 
method of active sets, originally developed for the discrete Cucker-Smale model, to the continuum de- 
scription. Their result on asymptotic fiocking can be easily adapted to our situation under the assumption 
that the communication rate function g is strictly positive: 



15 



Proposition 2 Let g{s) > go > for all s G [0, 1]. Consider the system (4.6) -(4.7) subject to compactly 

supported inital data {p^, u^) in M.'^ x M"^ and assume that it admits global smooth and compactly supported 
solutions {p{t),u{t)). Then the position and, resp., velocity diameters of the solution {p{t),u{t)), 

dx{t) = swp{\x - y\, x,y e supp p{t)}, 

du{t) = swp{\u{t,x) - u{t,y)\, x,y e supp p{t)}, 

satisfy 

sup dx{t) < +00, and d„(t) < du{0)e~^o^. 
t>o 

Proof: A slight modification of the proof of Proposition 6.4 in [17]. 

■ 

As already mentioned in Section 2, the assumption that g{s) > go > excludes the interesting case 
when there is no direct communication between certain regions of supp p, i.e., g{a[p]{x, y)) = for x G Jli, 
y & ^2, with some disjoint f^i, 0,2 C supp p. This immediately leads to the question if it is possible to 
modify the method developed in Section 3 to prove asymptotic flocking for the hydrodynamic system 
(4.6)-(4.7). The natural way would be to approximate the solutions {p,u) by finite systems of particles 
moving along characteristics, apply Theorem 1 for them, and pass to the limit, using appropriate stability 
properties. However, this program would be based on the a-priori assumption that the solution (p, u) is 
uniformly approximable by strongly connected systems of particles (in the sense of Theorem 1), which is 
a very strong and technical assumption. Moreover, the rigorous passage to the continuum limit for (2.1)- 
(2.2) would have to be developed, which seems to be, as mentioned before, far from trivial. Therefore, we 
will leave this topic open. 



5 Extension to other models of collective dynamics 

The idea of introducing the discreet (1.5) and continuum (4.3) topological distance is of course applicable 
to the full spectrum of models of collective behavior. Of course, this needs to be well justified from the 
modeling point of view, as in many cases it may be more appropriate to stay with the metric interactions. 

Let us provide just one particular example, the well-known attraction-repulsion particle model pro- 
posed in [10], 

Xi = Vi, (5.1) 
Vi = (a-6|i;i|2)^;,--i^VC/(|xi-x,|), i = l,...,N, (5.2) 

where a and b are nonnegative parameters and : M'^ — >■ M is a given potential modeling the short- 
range repulsion and long-range attraction. The term corresponding to a models the self-propulsion of 
individuals, whereas the term corresponding to b is the friction assumed to follow Raylcigh's law. The 
balance of these two terms imposes an asymptotic speed to the agent (if other effects are ignored), but 
does not infiuence the orientation vector. A typical choice for U is the Morse potential, which is radial 
and given by 

U{x) = -CAe-l^l/^^ + CRe-l^l/^«, (5.3) 



16 



where Ca, Cr and Ia, are the strengths and the typical lengths of attraction and repulsion, respectively. 

Clearly, the pairwise forces of attraction and repulsion decay with the metric distance between the 
respective agents. Let us now review the modeling assumptions of these two effects in the context of 
biological collective behavior: The repulsion force aims to describe the preference of the individuals to 
avoid collisions or even close contact (for instance in case of canibalism) with conspecifics. Therefore, it 
is a short-range metric interaction. On the other hand, the attraction force aims to model the preference 
of the individuals to form a compact group, e.g., a school, flock or swarm. The typical motivation for this 
behavior is to avoid predation [22]. Therefore, the strength of the attractive force should not be primarily 
related to the actual physical distance between individuals, and it seems appropriate to introduce the 
topological distance to model the attraction. For this we split the interaction term in (5.2) into two parts: 
the repulsion, still modeled as a gradient of the repulsion potential 

[/R(ar) = CRe-N/^«, (5.4) 

and the attraction, which, however, is not a potential gradient any more. Instead, we propose the new 
form 




i = l,...,N, (5.5) 



where qa ■ [0, 1] [0, oo) describes the strength of the attractive force in dependence on the topological 
distance a{xi,Xj) given by (1.5). The formal kinetic limit of (5.1), (5.5) reads then 

^ + ^; ■ V^/ + • [(a - b\v\^)vf] - • [{V^Ur * p)f] - V„ • {g[f]f) = 0, 
with p{t,x) = Jj^d f{t,x,v) dv and 

Q[m,x)= [ g{a[f]ix,y))^^p{t,y)dy, 
jRd \x - y\ 

with a[f]{x,y) given by (4.3). 

It is well known that the self-propelled attraction-repulsion model (5.1)-(5.2) produces a variety of 
different patterns, called clumps, ring clumping, rings and mills in [10]. It is out of scope if this paper to 
thoroughly investigate the patterning behavior of the modified model (5.1), (5.5). Let us only mention 
that our preliminary numerical simulations suggest that the new model is not only capable of reproducing 
all the above mentioned patterns, but also exhibits interesting novel patterning phenomena. In particular, 
we ran a simulation with N = 100 agents, initially randomly placed into the box (0, 1) x (0, 1) C M? with 
zero initial velocities. We used the repulsion potential (5.4) with Cr = 1.0 and £r = 0.5. The attraction 
in (5.5) was modeled using gA{s) = with Ca = 1-0 and Ia = 10/N = 0.1; this means that on 

average, each agent is mostly attracted by its 10 closest neighbors. Note that this scaling ensures that 
the repulsive and attractive forces are initially, i.e., when all the agents are randomly distributed in the 
unit box, of the same order of magnitude. We ran the simulation in the time interval t G [0, 10] and 
recorded the trajectories of all the agents in the upper left panel of Fig. 5. The agents quickly formed 
four groups of approximately same sizes and after an initial unordered movement, the groups started to 
follow an approximately circular path, being imposed by the interaction between the groups. However, on 
the scale of interactions within each group, we observe helixoidal-like movement of each single agent, as 
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Figure 5: Simulation of the model (5.1), (5.5) with N = 100 agents subject to a random initial condition 
in positions and zero initial velocities. Upper left panel shows the trajectories of all agents during the 
time interval t £ [0, 10]. A single agent trajectory for t G [8.0,8.5] is shown in the upper right panel, with 
another close-up in the lower left panel. Finally, another run of the simulation with a new random initial 
condition during the time interval t S [0, 100] generated the pattern shown in lower right panel, where 
again trajectories of all agents are recorded. 

depicted in the upper right panel of Fig. 5 for the time interval t € [8.0, 8.5] and another close-up in the 
lower left panel. It seems that this two-scale dynamic is genuine to the model with topological attractive 
interactions and as far as we know, cannot be produced by the original model (5.1)-(5.2). Interestingly, 
the quasi-circular pattern on the scale of the groups may not be persistent and may collapse into a chaotic 
movement, such as the one shown in the lower right panel of Fig. 5. This last figure was produced by 
running the simulation with all the same parameters, starting from another random sample for the initial 
condition, again with zero initial velocities. We conclude that the modified model (5.1), (5.5) exhibits 
very interesting and new patterning dynamics, which deserve further investigations. This is however not 
an objective of the present paper. 
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